Ferroelectric ordering in chiral smectic C* liquid crystals 
determined by nonchiral intermolecular interactions 
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General microscopic mechanism of ferroelectric ordering in chiral smectic C* liquid crystals is 
considered. It is shown that if the mesogenic molecules have a sufficiently low symmetry, the 
spontaneous polarization is proportional to one of the biaxial vector order parameters of the 
smectic C phase. This order parameter may be determined by intermolecular interactions which 
are not sensitive to molecular chirality. At the same time, the polarization is also proportional to 
a pseudoscalar parameter which vanishes if the molecules are nonchiral. The general statistical 
theory of ferroelectric ordering is illustrated by two particular models. The first model is 
based on electrostatic quadrupole-quadrupole interactions, and it enables one to obtain explicit 
analytical expressions for the spontaneous polarization. In the second model, the molecular 
chirality and polarity are determined by a pair of off-center nonparallel dipoles. For this case, the 
spontaneous polarization is calculated numerically as a function of temperature. The theory pro- 
vides a more general interpretation of the previous approaches including the classical Boulder model. 
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I. INTRODUCTION 

Ferroelectric smectic liquid crystals are unique sys- 
tems where the spontaneous polarization is determined 
by molecular chirality. Ferroelectric ordering in the chiral 
smectic C* phase has been predicted theoretically about 
three decades ago and both ferro- and antiferroelec- 
tric phases continue to attract a significant attention be- 
cause of their unusual structure and properties, and be- 
cause of their applications in electro-optical devices [2]- 
Ferroelectricity is observed only in tilted smectic phases. 
In these systems the spontaneous polarization is induced 
by the tilt and does not appear self-consistently like in 
conventional solid ferroelectric materials. In every tilted 
layer of a chiral smectic phase, the polarization appears 
in the direction of the polar C2 symmetry axis which is 
perpendicular to the tilt plane. In the bulk tilted phase 
molecular chirality also results in the formation of the 
macroscopic helical structure. In this structure, the di- 
rection of the tilt rotates while moving along the z-axis 
which is perpendicular to the smectic layers. As a result, 
chiral tiled smectics are characterized by the helical dis- 
tribution of the spontaneous polarizations and thus may 
also be called 'helielectric'. One notes that recently fero- 
and antiferroelectric ordering has been found in a novel 
class of smectic liquid crystal phases formed by achiral 
bent-core molecules. In those phases the spontaneous po- 
larization is not induced by the tilt and thus appears also 
in orthogonal smectic phases ^ . 

In the synclinic smectic C* phase, the direction of the 
tilt in adjacent layers is practically the same and thus the 
spontaneous polarization only slowly varies from layer to 
layer. In contrast, in the anticlinic smectic phase 
the direction of the tilt alternates from layer to layer to- 
gether with the polarization creating the structure with 
an antiferroelectric type ordering. In addition, many chi- 
ral smectic materials exhibit a sequence of the so-called 



intermediate phases in a narrow temperature interval be- 
tween the synclinic ferroelectric Smectic C* and the an- 
ticlinic antiferroelectric Smectic phase. Intermediate 
smectic phases are characterized by a 3D chiral distri- 
bution of the spontaneous polarization within the unit 
cell of 3 or 4 smectic layers [21 • Recently, the remarkably 
wide intermediate phases have also been discovered in 
mixtures of synclinic and anticlinic smectics where they 
can exist in a broad temperature range of up to 50° 
Finally, ferroelectric ordering also exists in smectic I and 
F phases which are characterized by some in-plane posi- 
tional or hexatic order. 

In spite of all diversity of tilted smectic phases with 
ferro-, antiferro- and ferrielectric ordering, the underly- 
ing mechanism is always related to the induction of the 
polarization by the tilt in individual chiral smectic lay- 
ers. Complex structures with a polarization distribution 
along the direction perpendicular to the layers appear 
due to interactions between the molecules in different lay- 
ers (see, for example, @, |^) which are generally weaker 
than intermolecular interactions within the same layer. 

The detailed microscopic mechanism of ferroelectric 
ordering in tilted smectics which is responsible for the 
induction of the polarization by the collective molecu- 
lar tilt has been the issue of debate during the past two 
decades. In particular, the role of molecular chirality has 
not been completely clarified. On the one hand, there is 
a general agreement that the spontaneous polarization in 
tilted smectics cannot exist without molecular chirality, 
i.e. at least a part of molecules must be chiral. On the 
other hand the role of chiral intermolecular interactions 
remains unclear. Some molecular models of the ferroelec- 
tric smectic C* phase are based on the assumption that 
the spontaneous polarization is directly determined by 
appropriate interactions between chiral molecules includ- 
ing, for example, interactions with a molecular chiral cen- 
ter 0, 0, [13, EH [13 ■ These models have been developed 
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FIG. 1: (Color online) a) Simple models for nonchiral 
molecules of the C2h symmetry used by Wulf and Goosens 
[20l.[2H. b) Corresponding models for chiral molecules where 
chirality is determined by the additional dipole perpendicular 
to the molecular plane. 




FIG. 2: (Color online) Schematic of biaxial molecule in the 
curved binding site as assumed in the Boulder model 



using an analogy with cholesteric liquid crystals where 
the helical twisting power is determined b y ch iral (albeit 
nonpolar) intermolecular interactions [l3l.Tl3.fl5|. Other 
models are based on a different microscopic mechanism 
[H, [13, [ll [H, [13, [lilwhich also requires molecular 
chirality, but, at the same time, takes into consideration 
only nonchiral intermolecular interactions. Some of these 
models are rather qualitative, but the underlying micro- 
scopic mechanism is essentially the same. 

The first extended description of this mechanism 
has been given by Goosens [l^, [2l[ who considered 
the electrostatic interaction between model molecular 
quadrupoles composed of two antiparallel dipoles which 
are perpendicular to the long molecular axis ( see 
Fig. [T]). Goosens has shown that the polarization in 
the smectic C* phase may be induced by the nonchiral 
quadrupole-quadrupole electrostatic interaction provided 
the molecules possess the additional dipole in the direc- 
tion perpendicular to the molecular plane. The latter 
dipole is responsible for the molecular chirality in this 
simple model, and the spontaneous polarization is pro- 
portional to the magnitude of the dipole and to the novel 
order parameter of the smectic C 20] phase which is re- 
lated to the relatively low symmetry of the molecule pre- 
sented in Fig. [TJ One notes, however, that the papers of 
Goosens are focused into one particular model, and do 
not contain any general theory of ferroelectric ordering 
or general expressions for the spontaneous polarization. 

A more general and a very successful model has been 
proposed by the Boulder group [H, In the Boul- 

der model the molecules of the zig-zag shape (see Fig. [2]) 
are ordering in the so-called binding cites which have the 
same point symmetry as the smectic C phase itself. Then 
transverse molecular dipoles are ordered in the particular 
direction perpendicular to the tilt plane simply because 
the zig-zag molecule fits into the binding cite of the same 
shape only for a particular direction of the transverse 
dipole. The Boulder model has been successfully used 
to describe and predict the value and sign of the spon- 
taneous polarization for a significant number of chiral 
smectic materials. This indicates that the corresponding 
mechanism of the ferroelectric ordering may be predomi- 
nant at least for conventional smectics C* . One notes also 



that the symmetry of a zig-zag molecule is exactly the 
same as that of the molecule with two equal antiparallel 
dipoles considered by Goosens. The interaction between 
the molecule and the binding cite , which is responsi- 
ble for the polar order, is also nonchiral in nature be- 
cause the binding cite itself is nonchiral. Boulder model 
emphasizes the steric mechanisms of the ordering, but 
the idea behind the model is much more universal. It is 
shown in this paper that in the context of a rather general 
molecular-statistical theory the concept of the 'binding 
cite' corresponds to the average one-particle mean-field 
potential which is created by all other molecules of the 
medium, and which reflects the symmetry of the smectic 
C phase. 

Terzis, Photinos, Samulski et.al. have developed a 
similar model [l^. This detailed model is based 
on a mean-field -like one-particle orientational potential 
for each molecular segment and, similar to the Boulder 
model, involves a summation over the molecular confor- 
mations. Using this model Terzis et.al. have obtained 
good quantitative results for the spontaneous polariza- 
tion of a number of ferroelectric smectics C* [iQj . 

Existing microscopic models for ferroelectric smectics 
C* have played an important role in the development of 
the theory of such materials. At the same time, from 
the point of view of theoretical physics, these models 
are too detailed in terms of a molecular structure and 
a particular choice of a coupling with the macroscopic 
environment. As a result, the models only indirectly ad- 
dress some of the more general physical problems related 
to the origin of ferroelectricity in tilted smectic phases 
including, for example, the description of the symmetry 
properties of relevant model interaction potentials and 
an interplay between the spontaneous polarization and 
the order parameters of the smectic C* phase. 

In this paper we develop a general mean-field 
molecular-statistical theory of ferroelectric ordering in 
the the smectic C* phase based on the general mech- 
anism described above, i.e., we consider the ferroelec- 
tric ordering in chiral smectics caused by nonchiral in- 
termolecular interactions. We obtain a simple expres- 
sion for the model interaction potential which may be 
responsible for ferroelectric ordering in tilted smectics, 



and interpret the general microscopic mechanism of the 
ordering mathematically using the concept of the pseu- 
dovector order parameter. The results of the theory will 
be used to obtain explicit expressions for the sponta- 
neous polarization in the smectic C* phase composed 
of biaxial molecules with quadrupole-quadrupolc inter- 
action potential. The spontaneous polarization together 
with other order parameters of the smectic C* phase will 
also be calculated numerically for another two interaction 
model; Finally, we will consider in detail the molecular 
origin of the ferroelectric ordering in novel smectic ma- 
terials which have recently been investigated by Lemieux 
et.al. [H, [13, m, These materials, which are used 
as chiral dopants in a nonchiral smectic C host, possess 
a molecular structure which is rather different from that 
of conventional ferroelectric smectics. In particular, the 
molecular chirality is mainly determined by the chiral 
distribution of permanent dipoles within the molecular 
structure. It is shown that despite a different nature of 
intermolecular interactions, the spontaneous polarization 
in smectics C doped with such molecules is still deter- 
mined by the same general mechanism. 

II. SPONTANEOUS POLARIZATION AND 
ORDER PARAMETERS OF THE SMECTIC C* 
PHASE 

A. Coupling between polarization, tilt and chirality 

It is well known since the work by R. Meyer [13] that 
from the purely macroscopic point of view the ferroelec- 
tric ordering in tilted smectics is determined by the linear 
coupling between the polarization and the tilt in a chiral 
medium. In terms of the Landau-de Gennes expansion 
the free energy of the ferroelectric smectic C* phase can 
be expressed as (see, for example, [T^): 

Fc = FA + F{e) + —P! + Cp(Ps-^), (1) 

XI- 

where Fa is the free energy of the smectic A phase, F{Q) 
is the expansion of the excess free energy of the smectic 
C phase in powers of the tilt angle Q and the last two 
terms describe the contribution which depends on the 
spontaneous polarization P^. Here w = (n • k)(k x n) 
is the so-called pseudovector tilt order parameter of the 
SmC phase where n is the director and k is the smectic 
layer normal as shown in Fig. [31 Minimization of the free 
energy (1) yields the well known result: 

Ps = CpW = Cp(n • k)(k X n), (2) 

which indicates that the spontaneous polarization is pro- 
portional to the pseudovector tilt order parameter and 
the coupling constant Cp. 

One notes that the polarization is a polar vector while 
the tilt order parameter w is a pseudovector with dif- 
ferent transformation properties (i.e, in contrast to the 
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FIG. 3: (Color online) Spontaneous polarization Ps and the 
pseudovector tilt order parameter w in the SmC* phase. 



polar vector P^ it does not change sign under space in- 
version). Thus the linear relationship (2) between polar- 
ization and the tilt is only possible if the coupling con- 
stant Cp is a pseudoscalar (which also changes sign under 
space inversion). Then the product of the pseudoscalar 
Cp and the pseudovector w makes the polar vector like 
polarization. Pseudoscalar quantities are nonzero only 
in a chiral medium. They are proportional to molecu- 
lar chirality and change sign when all chiral molecules 
reverse their handedness. Thus one arrives at a well es- 
tablished conclusion that the spontaneous polarization in 
tilted smectic phases occurs only if at least a fraction of 
molecules are chiral. 

At the same time, the pseudovector order parameter w 
is nonzero also in the nonchiral smectic C phase. Indeed, 
w is invariant under all symmetry transformations of the 
smectic C phase including the reflection with respect to 
the tilt plane, which is a symmetry plane. This is re- 
lated to the transformation properties of a pseudovector 
different from those of a polar vector. In this case, w is 
invariant under a reflection with respect to the tilt plane 
because both vectors n and k are in the tilt plane and 
thus are not effected by the reflection. In contrast, the 
spontaneous polarization Ps , of course, changes sign un- 
der a reflection with respect to the tilt plane. One notes 
that this does not violate the linear relationship (2) be- 
cause the pseudoscalar parameter Cp also changes sign 
under the reflection. In a nonchiral smectic C phase the 
coupling constant Cp vanishes identically and the sponta- 
neous polarization does not appear. 



B. Microscopic interpretation 

The purpose of any molecular theory of ferroelectric 
ordering in tilted smectics is to establish a relationship 
between the general macroscopic description presented 
in the previous subsection and the molecular ordering 
on the microscopic level. An intuitive interpretation of 
the ferroelectric ordering in the chiral smectic C* phase 
[H, [H [H can be illustrated using Fig. [H 

For illustration only, let us consider a simple model of 
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FIG. 4: (Color online) a) Molecular orientation A and B, 
which correspond to the opposite directions of the molecular 
transverse dipole, are equivalent in the SmA phase, b) In the 
SmC* phase molecular orientation A is more energetically 
favorable than the orientation B. 



a chiral biaxial molecule represented as a rigid rod with 
two 'lateral groups' and a permanent dipole perpendicu- 
lar to the molecular plane as shown in Fig. 2^. Note that 
the lateral groups make the molecule biaxial while the 
chirality is determined by the transverse dipole. With- 
out this dipole, the molecule possesses a mirror plane 
and thus it is nonchiral. Now let us assume that the lat- 
eral groups have a tendency to point in the direction of 
the region between two adjacent smectic layers. One can 
readily see that in the smectic A phase (i.e. without any 
tilt) the two orientations of such a molecule, which cor- 
respond to opposite directions of the transverse dipole fi 
are energetically equivalent. Thus the macroscopic po- 
larization in the untilted smectic phase should vanish. 
In contrast, in the tilted phase the balance between two 
opposite directions of the transverse molecular dipole is 
violated because the the molecular orientation A is more 
favorable than the orientation B (see Fig. 2b). As a re- 
sult, the average molecular dipole does not vanish, and 
a macroscopic polarization appears in the direction per- 
pendicular to the tilt plane. Now one has to clarify how 
this type of ordering corresponds to the general macro- 
scopic description presented above. 

Firstly one notes that this simple argument is valid for 
molecules which are characterized by the C2h symmetry 
when the transverse dipole is removed. Simple examples 
of molecules of the C2h symmetry, presented in FiglU in- 
clude a molecule with two in-plane antiparallel dipoles 
(a model considered by Goosens [1^) and a molecule of 
the 'zig-zag' shape (considered by Wulf [1^ and others) 
which very roughly characterize the actual shape of typi- 
cal mesogenic molecules. Orientation of any rigid biaxial 
molecule can be specified by the unit vectors a and b 
in the direction of short and long molecular axis, respec- 
tively. In practice, the orientation of a nonpolar molecule 
is usually characterized by second rank tensors composed 
of the components of the vectors a and b. For exam- 
ple, uniaxial molecules are characterized by the molecular 
tensor — 0^0^ — (l/3)(5„/3- The statistical average 
of Q^^ is the nematic tensor order parameter Q. Biaxial 
molecules of high symmetry, which possess two mutually 
perpendicular symmetry planes, are also characterized by 
the second molecular tensor = habp — CaCp where the 



unit vector c_Lb is the second molecular short axis. In 
addition, the molecules of the C2h or lower symmetry are 
characterized by the third molecular tensor 0^6^ which 
is invariant under all symmetry operations which leave 
the molecule intact. Indeed, one can readily see that the 
molecules presented in Figs (T] and [5] are not invariant un- 
der the sign inversion of the axes a or b individually. At 
the same time, the molecules are invariant under simul- 
taneous inversion of both axis a and b. This symmetry 
enables one to introduce the transverse molecular pseu- 
dovector (ax b) which is related to the skew part of aQ,6^. 
One notes that the existence of this transverse pseudovec- 
tor does not violate the mirror symmetry of the molecule 
because the pseudovector (a x b) is invariant under a re- 
flection with respect to the molecular mirror plane which 
is parallel to a and b. 

Now it can readily be shown that the average (a x b) is 
nonzero only in a tilted smectic phase and is proportional 
to w = (n • k)(k X n) Thus the expression 



v^ab = (a X b) 



(3) 



is a microscopic definition of a pseudovector tilt order 
parameter of the smectic C phase. One notes that in 
the smectic C phase composed of biaxial molecules there 
exist several tilt order parameters. However, only the 
parameter ([3]) is directly related to the spontaneous po- 
larization. 

Indeed, one notes that the molecular orientations A 
and B in Fig. [5] are characterized by the opposite direc- 
tions of molecular pseudovector (a x b) (because the the 
short axis b points in the opposite directions while the 
long axis a is the same). Thus one concludes that the 
average pseudovector (a x b) 7^ in the smectic C phase 
where the orientation A is more energetically favorable 
than the orientation B. Moreover, it follows from the gen- 
eral symmetry arguments that (axb) c>c w = (n-k)(kxn) 
because w is the only pseudovector allowed by the sym- 
metry of the smectic C* phase. Indeed, any macroscopic 
vector or pseudovector must be parallel to the C2 sym- 
metry axis of the smectic C* layer which is normal to the 
tilt plane. This is exactly the direction of w. A different 
derivation of this result is presented in Section 3. 

Finally it can be shown that the spontaneous polar- 
ization Ps is proportional to the average (a x b) oc w. 
We assume for simplicity that the transverse molecular 
dipole is parallel to the short molecular axis c_Lb_La. i.e 
jj, = fj,±c. Then the macroscopic polarization in the 
smectic C* phase equals 



(4) 



where p is the molecular number density. 

The unit vector c can be expressed in terms of the unit 
vectors a and b in the following way 



c = A(a X b) 



(5) 



where A — ((ax b) -c) is the molecular unit pseudoscalar 
which specifies the handedness of the molecular coordi- 
nate system. Note that in Eq.(5) c is the conventional 
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polar vector which is expressed as a product of the pseu- 
dovector (a x b) and the pseudoscalar A. 

Now the spontaneous polarization (4) can be rewritten 
as: 

Ps = PM-l(c) = PM±A(a X b) = pA^(a x b), (6) 

where = (/x-(axb)). Here we have taken into account 
that the parameter A is independent of the molecular 
orientation. 

According to Eq.(6), the spontaneous polarization 
is proportional to the tilt order parameter (a x b) cx w = 
(n • k) (k X n) . Thus the simple microscopic interpreta- 
tion of the appearance of the spontaneous polarization 
in the smectic C* phase, presented in this subsection, is 
fully consistent with the general phenomenological the- 
ory outlined in subsection III Al 

In summary, one concludes that if the smectic C* phase 
is composed of molecules of the C2h symmetry or lower 
with an additional transverse dipole in the direction per- 
pendicular to the molecular plane, the spontaneous polar- 
ization is proportional to the order parameter ^ and the 
pseudoscalar quantity A^ = {/i- (ax b)) which character- 
izes the molecular chirality determined by the transverse 
dipole fi. A consistent statistical theory of the ferroelec- 
tric ordering in the smectic C* phase is presented in the 
following section. 



III. MOLECULAR STATISTICAL THEORY OF 
FERROELECTRIC ORDERING IN THE 
SMECTIC C* PHASE 

A. General results 

As discussed in the previous section, the orientation of 
a biaxial molecule can be specified by the unit vectors a 
and b in the direction of the long and short molecular 
axis, respectively. The second molecular short axis c is 
then given by Eq.(5). The spontaneous polarization can 
be expressed in the following general form: 



J /x_L/i(a, b)dadb. 



(7) 



where /i (a, b) is the orientational distribution function of 
the smectic C* phase, and the transverse molecular dipole 
fj,± = /i^c. The orientational distribution function can 
always be expressed in terms of the effective one-particle 
potential f/i(a, b): 



/i(a,b) 



1 



■ exp 



f/i(a,b) 



where 



Z 



cxp 



t/i(a,b) 



da dh. 



(8) 



(9) 



For molecules which are polar in the direction of the 
c-axis, in the quadrupolar approximation the effective 



one-particle potential depends on the unit vector c and 
the second rank molecular tensors a^ap, bahp and ajjp 
which are invariant under all symmetry operations of a 
molecule. Here c is expressed in terms of a and b by 
Eq.(5). As discussed in subsection IH Bl the invariant 
aabp exists only for biaxial molecules of sufficiently low 
symmetry. Thus in the quadrupolar approximation the 
potential [/i(a, b) reads: 

C^i(a, b) = a^apAap + babpBai3+ 

{aj)f3 + baap)Ca[3 + (a X b) • W, (10) 

where W and Aq,^, Baf3, Cap are the material pseudovec- 
tor and tensors correspondingly, which depend on the 
symmetry of the chiral smectic C* phase. This sym- 
metry, in turn, is determined by the two tensors UaUp 
and kakp which specify the macroscopic structure of the 
phase. One notes that the last term in Eq.(lO) contains 
a pseudovector (ax b) , which means that W has also to 
be a pseudovector in order to insure that C/i(a, b) is a 
scalar. We conclude that this last term is determined by 
chiral interactions in the system, i.e., the interactions de- 
pending on the handedness of the interacting molecules. 
Indeed, let us consider the case when Eq.(lO) describes 
an effective one-particle potential of a chiral dopant in 
the chiral smectic C* host. Inversion of the host chirality 
results in sign inversion of the pseudovector W while all 
material tensors Aap, Bap and Cap remain the same. As 
a result, the last term in Eq.(lO) changes sign. Thus the 
last term describes the so called chiral discrimination en- 
ergy, i.e. the difference of energy of interaction between 
the same chiral molecule and the two enantiomeric forms 
of the chiral host. In many cases such a discrimination 
is small and then this term may be neglected. 

Taking into account that all macroscopic quantities in 
the smectic C phase must be quadratic both in n and 
k, the pseudovector W can be expressed in the following 
general form: 



W 



KW, 



(11) 



where w is a pseudoscalar and k is a scalar. The sym- 
metric tensors A, B and C are expressed as: 



AaP = AiHaUp + A2kakp + ^3(11 • k)(n„fc^ -t- kaUp), 



Bap = BmaUp + B2kakp -|- ^3(11 • k){nakp + kaup), 



Cap = Cinanp+C2kakp+C3{n ■ k){nakp+kanp), (12) 

where we have neglected the terms quadratic in and 
w because these terms are of the order of at small tilt 
angle 8^1 while all other terms in Eq.(12) are of the 
order of 1. 

The effective one particle potential Ui (a, b) is deter- 
mined by intermolecular interactions in the smectic C* 
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phase which may, or may not be sensitive to molecu- 
lar chirality. For example, the electrostatic interaction 
between permanent dipoles and quadrupoles is not sen- 
sitive to molecular chirality while the interactions in- 
volving molecular octupoles arc different for the pairs 
of molecules of equal and opposite handedness, respec- 
tively. In particular, the pseudovector quantity W in 
Eq.(ll) must vanish if the molecules are nonchiral. This 
parameter is determined by some chiral intermolecular 
interactions which exist only between chiral molecules. 
All other material parameters in Eqs. (10-12) are scalars, 
and therefore they are generally nonzero in the corre- 
sponding nonchiral smectic C phase. The difference be- 
tween scalar and pseudoscalar material parameters in the 
effective one-particle potential enables one to distinguish 
between two different microscopic mechanisms of ferro- 
electric ordering in the smectic C* phase. 

Taking into account that the spontaneous polarization 
Ps and the pseudovector w are small at small tilt angles 
Q, the orientational distribution function (8) can be ex- 
panded in powers of and w keeping the linear terms, 
and substituted into the general expression for the spon- 
taneous polarization (7). As a result, the spontaneous 
polarization in the smectic C* phase can be expressed 
as: 



^ • Ps = pAnX ■ w- 



1 

/^i^exp 

^0 



^r(a,b) 



da dh, (13) 



where the inverse polarizability tensor = 
Saf3 + {p/kBT){u)^i±aCf3)o, the tensor Xap = 
~ip/kBT){p±aC(})o , the averaging (...)o is per- 
formed with the orientational distribution function 
/o = (l/Zo)exp(-C/}°)(a,b)/fcBr) and C/l°^(a,b) is 
given by Eqs. (10-12) with W = 0. 

According to Eq.(13) there exist two qualitatively dif- 
ferent contributions to the spontaneous polarization of 
the smectic C* phase which correspond to the two terms 
in the r.h.s. of Eq.(13). The first term comes from the 
last term in the one-particle potential (|10p . which is de- 
termined by chiral intermolecular interactions. This con- 
tribution corresponds to the microscopic mechanism of 
ferroelectric ordering determined by chiral interactions 
which has been considered in refs.fg, [l^, [HI, [I^l- I^i con- 
trast, the second term describes the ferroelectric ordering 
of chiral molecules determined by nonchiral intermolecu- 
lar interactions. Indeed, this contribution describes the 
ordering of a transverse molecular dipole in a nonchiral 

effective one-particle potential u[^\a,h). This potential 
does not depend on any pseudoscalar parameters, and 
thus, in the first approximation, it is exactly the same as 
in the corresponding nonchiral smectic C phase. Molecu- 
lar chirality in this case is determined by the orientation 
of the molecular transverse dipole and manifests itself 
during the averaging process as described below. Now 
let us consider this contribution in more detail. 



In the general case the effective one-particle potential 
u[^^ (a, b) can be expressed as a sum of the following two 
terms: 



C/r(a,b) 



Uaia) + Ubia, b). 



(14) 



where Ua{a) — aaa/^Aap and Ub{a,h) = hahpEap + 
[o-abfj + baap)Cai3- In Eq.(14) the first term depends 
only on the orientation of the long molecular axis a while 
the second term depends also on the orientation of the 
short axis b. Smectic liquid crystals are composed of 
strongly anisotropic molecules, and in this case one may 
assume that the intermolecular interaction energy associ- 
ated with a change of orientation of the long axes is much 
larger than the change of the energy associated with the 
rotation of short molecular axes. Then the second term in 
Eq.(14) is expected to be much smaller than the first one, 
i.e. Ub <C Ua- One notes that the first term in Eq.(14) 
determines the orientational (nematic) ordering of long 
molecular axis, while the second term is responsible for a 
weaker ordering of short molecular axes. Now the orien- 
tational distribution function /o(a, b) can be expanded 
in powers of Ub keeping the linear term: 



1 

:^exp 



kuT 



1 - 



Ub{a, b) 



knT 



, (15) 



where the uniaxial orientational distribution function 
/a (a) depends only on the orientation of the long axis 



/a(a) 



1 



cxp 



Ugja) 
keT 



(16) 



Substituting Eq.(15) into the second term in Eq.(13) one 
obtains the following expression for the spontaneous po- 
larization 



Ps ~ A 



'ksT 



J {a X b) Ubia,h) faia)dadh. (17) 



Now Eqs. (14) and (16) can be substituted into Eq.(17) 
where the averaging is performed over b and a taking 
into account that the function /a (a) is independent of b. 
Neglecting biaxiality of the smectic C phase one may use 
the following simple formulae: 

If 1 

— babjsdh = -{5ai3 - aaap); (18) 

and 

If 1 1 

— / fa{a){aaai3 - -Sap) da = S{nanp - -Sap)-, (19) 

and obtain the final expression for the spontaneous po- 
larization: 



n X k](n • k). 



2kBT 

where S is the nematic order parameter. 



(20) 
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Eq.(19), which has been obtained without using any 
particular model, presents a general expression for 
the spontaneous polarization of the chiral smectic C* 
phase determined by nonchiral intermolecular interac- 
tions. The spontaneous polarization is proportional to 
the pseudoscalar parameter = (/^±(a x b)) which 
specifies molecular chirality determined by the orienta- 
tion of the transverse dipole, and the parameter C3 which 
is determined by nonchiral interactions. One notes that 
the general microscopic mechanism of ferroelectric order- 
ing described by the present theory is qualitatively the 
same as the one considered in the Boulder model [1^, [13] ■ 
In the framework of the Boulder model the restrictions 
imposed by the binding cite on the orientation of a given 
molecule are equivalent to the specific form of the effec- 
tive one-particle potential t/}"'' (a, b) . From the qualita- 
tive point of view, both the binding cite and the effective 
one-particle potential possess the same symmetry as the 
nonchiral smectic C phase (i.e. the C2 symmetry axis 
and the mirror plane) and are not sensitive to molecular 
chirality. On the other hand, the present theory also cov- 
ers the model of Goosens .21] based on the electrostatic 
interaction between molecules with pairs of antiparallel 
dipoles. Detailed numerical results obtained using a ver- 
sion of this model are presented in Section IIVI Explicit 
expressions for effective one-particle potential and the 
spontaneous polarization in the quadrupole-quadrupole 
interaction model are derived in the following section us- 
ing the mean-field approximation. 



B. Quadrupole-Quadrupole interaction model 

In this section we obtain explicit analytical results for 
the spontaneous polarization in the smectic C* phase us- 
ing the simple model of a rigid molecule with essentially 
uniaxial quadrupole tensor slightly tilted with respect to 
the primary molecular axis (which is determined by the 
shape or molecular inertia tensor and orders along the 
director). Such a molecule is biaxial, and the molecular 
biaxiality is determined by the angle between the axis of 
the quadrupole tensor and the primary molecular axis. 
Assuming that this angle is small, one obtains approxi- 
mate expressions for the polarization using the general- 
ized mean-field theory. 

Let us consider the system of rigid molecules with pri- 
mary axes SLi. The molecules are characterized by the 
permanent quadrupole tensor which depends on the 
distribution of effective charges within a molecule. It 
should be noted that for molecules of the C2h symme- 
try, considered here, one primary axis of any molecular 
tensor (including the quadrupole one) must be parallel 
to the two-fold symmetry axis of the molecule c. At the 
same time, the orientation of the two remaining primary 
axes, which lie in the symmetry plane of the molecule, is 
not specified by the molecular symmetry. As a result, the 
orientation of these axes should generally be different for 
different molecular tensors characterizing different molec- 



ular properties. In particular, the primary axes of the 
quadrupole tensor Qq,^ are not expected to coincide with 
those of the molecular inertia tensor, a and b. Thus in 
the general case the traceless symmetric quadrupole ten- 
sor can be expressed in terms of molecular axes a_Lb_Lc 
in the following way: 

QaP = qiiaaap-Sal3/i)+q2{babp~CaCp)+q3{aabi3 + baaf3). 

(21) 

It will be shown below that the spontaneous polarization 
is proportional to the off-diagonal element 93, which char- 
acterizes the difference in the orientation of the molecular 
inertia and quadrupole tensors. 

The electrostatic interaction between permanent 
quadrupoles of the molecules i and j can be written in 
the form: 

UqqiiJ) = QapDapysqi^s: (22) 

where the quadrupole-quadrupole coupling tensor D is 
given by 

3 

5Saf3UjUS — bSajUpUg — 5SaSU-yUf3 — 56f3jUsUa — 

56si3U-yUa — 5S-ySUaUf3 + 35UaUpU^Us), (23) 

where R = Ry is the intermolecular vector and the unit 
vector u = R/|R|. 

Now the quadrupole-quadrupole interaction can be 
taken into account in the generalized mean-filed approx- 
imation of the smectic C* phase. In this approximation 
(see, for example [s^, HH [13] ) the free energy of the sys- 
tem without positional order can be expressed as: 

pksT J hini)lnfi{ni) dn^, (24) 

where /i(r2) is the one-particle orientational distribution 
function, the variable Qi — [ai,bi] specifies the orienta- 
tion of the molecule i, U{1, 2) is the effective pair interac- 
tion potential which takes into account the steric cut-off, 
and p is the number density of molecules per unit area 
of the smectic layer. 

Minimization of the free energy (24) yields the one- 
particle distribution fimction in the form ([8l [9]) with the 
one-particle potential Ui being equal to the mean- field 
potential 

UmfM =pJ C/(l,2)/i(a2,b2)d2R dh2 dsi2- (25) 

The total effective pair interaction potential U{1,2) 
can be expressed as 

;7(l,2) = ?7o(l,2) + [/,,(l,2), (26) 

where C/o(l,2) is the effective interaction potential for 
molecules without permanent quadrupole moments, and 
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Uqq{l^ 2) is the quadrupole-quadrupole interaction energy 



given by Eq. (|^^ . We assume for simplicity that the po- 
tential Uo{l, 2) is even in ai, a2, bi, b2,, and thus it can- 
not be responsible for the ferroelectric ordering in the 
smectic C* phase. 

Substituting Eqs. (l22l23|) for the quadrupole- 
quadrupole interaction potential into Eq. ()26p and 
then into Eq. pS)) one obtains the following expression 
for the total mean-field potential after averaging over 
a.2 , b2 and integration over R: 



UMF{a, b) = C/|?^(a, b) -I- Ul^pia, b), 
where f7|^^(a, b) is even in a and b, and where 



(27) 



(28) 



Here q^^l is the quadrupole tensor of the molecule 1 given 



by Eq. ()2ip . Rq is the distance of minimum approach be- 
tween the neighboring molecules in the layer, and the 
tensor D is given by Ea. ((23| with R — Rq and u = k. 

At this point we can neglect the weak biaxiality of 
the molecular distribution and use the average of the 
quadrupole tensor pTjl in the simple form 



(29) 



(gap) = qiS {Uanp - Sap/S), 



After some algebra the corresponding expression for 
DofySp{'i<-){q'ys) can be written as: 



the biaxial ordering in the system is weak (i.e. using 
Eqs. (|18ll9p for the averaging) one obtains the following 
expression for the spontaneous polarization of the smec- 
tic C* phase: 



pA,,c;;[nx k](n-k), 



with 



ISvrp 



rnksTRl 



(?ig3^'(4-7sin2e). 



(34) 



(35) 



This explicit analytical expression for the spontaneous 
polarization was obtained in the context of the model 
of the quadrupole-quadrupole intermolecular interactions 
and confirms all results of the previous section which have 
been obtained using general theory arguments. In partic- 
ular, the spontaneous polarization is proportional to the 
pseudoscalar parameter = (/x • [a x b]) which speci- 
fies molecular chirality. At the same time, the remaining 
factor Cp is completely independent of the molecular chi- 
rality, and is obtained using the orientational distribution 
function of the nonchiral smectic C phase. In this simple 
model, the factor depends only on the diagonal com- 
ponents of the molecular quadrupole moment, the orien- 
tational order parameters of the smectic C phase and the 
distance of minimum approach Rq. 

In the following section we consider a more realistic 
model based on electrostatic interactions between local- 
ized molecular dipoles. 



Ma/3 = Da(3'ysO<.){q-fS) 



[Sapi^ — 5cos^ 6) + 2nanp — 5(1 — Tcos^ Q)kakp— 

10(k-n)(n„fc/3 + n/5fca)]. (30) 

Assuming that the off-diagonal element 93 of the 
molecular quadrupole tensor is small one may expand the 
exponent in Eqs. (j8l9p in powers of 53 keeping the linear 
terms. This yields the following approximate expression 
for the orientational distribution function: 



/i(a,b) -/o(a,b) + A/(a,b), 



(31) 



where 



/o(a,b) = — exp 
^0 



C/i°^(a,b) T:pRlq%Map 



ksT 



AkeT 



and 
A/(a,b) 



(32) 



■Kp 

^keT 



Rlq-i{aahp + apha)Map fo{ai,h), (33) 



and where is the molecular quadrupole tensor given 
by Ea. ([2T|) without the off-diagonal term. 

Only the term containing A/ makes a contribution to 
the spontaneous polarization given by the general equa- 
tions (|6l7p . Substituting Eas. (|31ll55|) and assuming that 



IV. INTERACTION BETWEEN CHIRAL PAIRS 
OF MOLECULAR DIPOLES 

In a recent series of papers by Lemieux et.al. [H, 
[25l . [2^ a number of novel compounds with unconven- 
tional structure have been used as chiral dopants to in- 
duce the large spontaneous polarization in the smectic C 
phase. One notes that the chirality of these molecules 
is not determined by any chiral centers, but is a con- 
sequence of a chiral distribution of permanent molecular 
dipoles. Disregarding other elements of the actual molec- 
ular structure, one can use the minimum model shown in 
Fig. [5l In this model, the molecule is presented by a rigid 
rod (with some dispersion interactions between rods) and 
a pair of off-center dipoles with large transverse compo- 
nents lying in orthogonal planes. Introducing the orthog- 
onal transverse unit vectors we write 



^ = Af(e^ sina ± acosa). 



(36) 



One can readily see that the molecule presented in Fig. [5] 
is chiral because it does not have any symmetry planes. 
The total dipole moment of the molecule is transverse 
/x^ = (/x+ -f /x~) and directed along the unit vector 
c= (e++e^)/\/2. Accordingly, another short molec- 
ular axis is to be defined as b = (e+ — e~ )/ -^2. The 
spontaneous polarization in such a system is given by 
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V 



FIG. 5: Schematic of molecular axes and dipoles 



the general eqs. (|4l6p where the pseudoscalar parameter 
is expressed as: 



V2 
ji sin a 



a • X /X ] = \/2A /i sin ( 



(37) 



This parameter quantitatively determines the molecular 
chirality in the context of the present model.. 

An interaction potential for a pair of such molecules is 
expressed as a sum of effective interaction potentials be- 
tween rigid uniaxial cores and the sum of all electrostatic 
dipole-dipole interactions: 

f/(l,2) = C/'^"(ai,R,a2) + [/^(l,2) (38) 

with depending also on the orientation of short molec- 
ular axes: 

where the both tensors D involved have the form 



Ai(ai, a2, R) = {r^Sij - Snrj) 



(40) 



with the distance between the interacting dipoles r being 
the function of molecular orientation and intermolccular 
distance: r — r(ai,a2,R). For D|| it reads r = R + 
i^a2 — vELi, while for it is r = R + i^a2 -I- i'sli. 

We assume that the tilting of the director in the smec- 
tic C phase is due to the long- axes potential U'^'^{1,2). 
We then employ the following model expression for the 
uniaxial potential which has been extensively used in the 
general theory of Sm^-SmC transition 



C/^^(ai , R, a2) « z;i (R) [(ai • u)^ + (a2 • u)^] 
+ V2{R){ai ■ f + V3{R){ai ■ a2)(ai • u)(a2 • u) 

-K«4(i?)(ai •u)2(a2 •u)2, (41) 

As discussed in detail in [H, the corresponding mean- 
field potential 

U^p{a) = wiP2(cos7) + W2SkP2{cosj)+ 

wzPk sin^ 7 cos 2(j) + WiV sin 27 cos 0, (42) 



depends on three order parameters 

Sk = {P2 (cos 7)) , Pk = (sin2 7 cos 2^) , (43) 
V = (sin 27 cos tp) , (44) 

where 7 and are the polar and azimuthal angles of the 
unit vector a respectively. 

If these order parameters are known, conventional or- 
der parameters such as nematic order parameter 5, ne- 
matic tensor biaxiality P, and the tilt angle can be 
easily calculated as 



tan 26 



V 



Sk - O.SPfc ' 
3 V 



1 3 

s = -^Sk + -Pk , 4^.^^2e' 

B = -Sk + -Pk - 2— 



(45) 
(46) 
(47) 



The potential ([55]) promotes the tilt if the nematic order 
parameter exceeds the critical value 



Sac 



3wi 



4u'4 — 3W2 ' 



(48) 



which means that the growth of the nematic order is 
the driving force of the tilting transition. Thus in this 
model the molecular dipoles are not responsible for the 
tilt of the director, but the interaction between such pairs 
of dipoles gives rise to the spontaneous polarization as 
shown below. 

One notes that the actual form of the interaction po- 
tential for off-center dipoles is too complicated to be 
used directly in the statistical theory |3J|. In particular, 
substituting the actual potential ([55]) into the Eg. ([^5]) 
one cannot obtain the mean-field potential as an ex- 
plicit function of the orientational order parameters. As 
a result the free energy cannot be minimized to deter- 
mine the transition point and the temperature variation 
of the parameters. However, we can expand the actual 
dipole-dipole potential in spherical invariants neglecting 
the higher order terms. One notes that the statistical 
averages of higher order terms are expressed in terms 
of higher order orientational order parameters which are 
normally not important from the qualitative point of view 
[3^ . The details of the expansion procedure are pre- 
sented in the Appendix \X\ 

As shown in the Appendix |^ the actual interaction 
potential between the pairs of off-center dipoles can be 
approximated by the relatively simple expression: 

C^MF = w^lcota V{a-k){h ■ t) + 

cot a r sin 27 cos </? -I- \/2 T (a • k)(b • t)], (49) 

which involves the biaxial order parameter T = ((a-k)(b- 
t)). Combining this with the uniaxial mean-field poten- 
tial (j42p we can write the total mean-filed approximation 
of the LC free energy ([M)) . 
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FIG. 6: (Color online) (a) Smyl-SmC phase transition ob- 
tained using the parameters wi = —0.05, W2 = —1, W3 = 
—0.9, W4 = —0.8, W5 = (thin lines), kjs = —0.05 (dashed), 
Wh = —0.1 (solid) and q = 30°. Biaxial order parameters V 
and Gk together with the normalized polarization Ps/ p/^f_i are 
shown in (b) for the case wz, — —0.05 and in (c) for ws — —0.1 



Now the orientational order parameters can be evalu- 
ated by numerical minimization of the free energy at a 
given temperature, and the spontaneous polarization can 
be calculated using the general Eq.(I7]) for the orienta- 
tional distribution function ([5]) and the sum of Eqs. (UH) 
and (|49)) as the total mean-field potential. Typical results 
of these calculations are presented in Fig. [51 

Thus the present general theory enables one to calcu- 
late the spontaneous polarization explicitly as a function 
of temperature using the model of rod- like molecules with 
chiral distribution of dipoles. In this case the ferroelec- 
tric ordering also occurs according to the general mecha- 
nism described in the previous sections. Indeed, one can 
readily see from Fig. [5] that the spontaneous polariza- 
tion is approximately proportional to the order parame- 
ter r which is mainly determined by the nonchiral part 
Uo{l,2) of the total interaction potential. At the same 
time, the spontaneous polarization is proportional to the 
pseudoscalar parameter oc (a • [fj,~^ x which, of 

course, vanishes if the molecules are nonchiral, for ex- 
ample, if both dipoles and the long axis of the molecule 
presented in Fig. [5] are in the same plane. 



V. CONCLUSIONS 

I this paper we have used the general statistical theory 
and two particular molecular models to demonstrate how 
the ferroelectric ordering of polar and chiral molecules 
in the Smectic C* phase may be determined by nonchi- 
ral intermolecular interactions. It has been shown that 
if the molecules are characterized by C2h symmetry or 
lower, the spontaneous polarization is given by the gen- 
eral expression Pg = pAp(a x b) (see Eq.(IB])), which 
is an exact result derived without any approximations. 
Here the parameter = {fi± ■ (a x b)) characterizes 
the molecular chirality determined by the relative orien- 
tation of the transverse molecular dipole and the molecu- 
lar plane containing the long and short molecular axes, a 
and b respectively. For rigid molecules of C2/1 symmetry 
(see FigH]), the spontaneous polarization is always pro- 
portional to the pseudovector order parameter (a x b) 
which is nonzero also in the nonchiral Smectic C phase. 
Thus, if specific chiral and polar intermolecular interac- 
tions are not important, the spontaneous polarization is 
essentially determined by the molecular chirality coeffi- 
cient A^ and the order parameter (a x b) emerging due 
to nonchiral molecular interactions in the corresponding 
nonchiral Smectic C phase. 

One notes that this general mechanism of ferroelec- 
tric ordering corresponds to the one considered in the 
Boulder model 16, 17]. In the Boulder model, a sin- 
gle molecule of sufficiently low symmetry is ordered in 
the binding cite which plays the role of the effective one- 
particle (mean-field) potential considered in the present 
statistical theory. Similar to the effective mean-field po- 
tential, the binding site itself has exactly the same sym- 
metry as the nonchiral smectic C phase (see Fig. i.e., 
possesses the mirror plane and the two-fold symmetry 
axis and thus is nonchiral. One notes also that, in the 
first approximation, the ordering of a molecule of the 
zig-zag shape in the binding site is not determined by its 
possible chirality. The corresponding nonchiral molecule 
with a similar overall shape will also order in the same 
binding site, although no polarization will be created in 
this case, of course. The ordering of such a nonchiral 
molecule in the binding site is described by the same 
nonchiral pseudovector order parameter w which is con- 
sidered in this paper. The same general mechanism cor- 
responds also to the one considered in [l^, [l^, for 
some particular cases. At the same time, the Boulder 
model as well as the models considered in [3, enable 
one to account for a selection of molecular conformations 
which fit the site or minimize some interaction potential. 
This effect is not taken into consideration in the present 
paper which deals with rigid molecules. 

In this paper, the general mechanism of ferroelectric 
ordering has been illustrated using two particular molec- 
ular models. The first model is based on the electrostatic 
interaction between anisotropic molecular quadrupoles. 
Here the molecular quadrupole tensor is assumed to be 
nondiagonal in the molecular frame determined by the 
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molecular axes a, b and c (i.e. one of the primary axes 
of the molecular quadrupole is tilted with respect to the 
long axis a) . The existence of the nonzero off-diagonal 
element (see Eg. ipT]) ') determines the C2h symmetry 
of the molecule. If this off-diagonal element is small, it 
is possible to expand the orientational distribution func- 
tion in powers of and obtain the explicit analytical ex- 
pression for the spontaneous polarization proportional to 
(73. This models shows that the electrostatic quadrupole- 
quadrupole interaction, which is not sensitive to molecu- 
lar chirality, may be responsible for the ferroelectric or- 
dering of chiral molecules in the Smectic C* phase. 

Finally we have considered a more realistic molecular 
model related to the materials recently synthesized by 
Lemieux et.al. ^2^, [23, H^, HH] , in which the molecular 
chirality is mainly determined by the distribution of per- 
manent dipoles. We consider a simple model of a uniaxial 
rod with two nonparallel off-center dipoles which make 
approximately an angle of 90°. In this simple case, the 
pair of permanent dipoles is responsible for both molec- 
ular chirality, polarity and biaxiality. The model inter- 
action potential for such molecules is composed of the 
uniaxial interaction responsible for the tilt in the smec- 
tic C phase, and the electrostatic interaction between all 
dipoles. We obtain that here the ferroelectric ordering 



also follows the general mechanism described above, and 
the spontaneous polarization is proportional to the pseu- 
dovector order parameter w. The polarization is also 
proportional to the pseudoscalar parameter which 
vanishes if we set the molecule to be nonchiral by placing 
the two dipoles and the long molecular axis a within the 
same plane. One notes that for real materials of the type 
reported in [23,, 24, 25, 26], the direct interaction between 
pairs of dipoles may not be the only cause of spontaneous 
polarization. Such molecules may also possess a zig-zag 
shape or have conformational states of the corresponding 
symmetry. Then the molecule would order in the bind- 
ing site according to the Boulder model just due to steric 
interactions, and this will make an additional contribu- 
tion to the spontaneous polarization. Which contribu- 
tion is predominant for particular materials can be de- 
termined by experiments involving systematic variation 
of the molecular structure. 
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APPENDIX A: APPROXIMATION FOR THE 
ELECTROSTATIC INTERACTION BETWEEN 
PAIRS OF DIPOLES 

Let us consider the electrostatic dipole-dipole interac- 
tion p9|) between two molecules described in Section lTVl . 
In order to calculate the corresponding part of the mean- 
field potential ([^5]) , one has first to integrate ([55)1 over the 
intermolecular distances. This involves the integrals 

D(ai,a2,k)=/ D(ai, aa, R) d^R. (Al) 

jR^n(ai,a2) 

taken over the distances R within the smectic plane and 
accounting for the steric cut-off between the rigid cores 
of the molecules with long axes ai and a.2. 

The tensors (|A1[) have important properties: they are 
invariant under the following transformations 

ai ^ a2 and ai, a2 ^ -ai, -a2, (A2) 

and are traceless and symmetric: 

(5y = 0, D,, = Dj,. (A3) 

Furthermore, the substitutions 

ai,a2^-ai,a2 and ai,a2^ai,-a2 (A4) 

transform into D|| and vice versa. 

We are interested in the orientational interaction which 
containing the lowest possible (i.e. first) power of ai_2- 
There are four tensorial expressions obeying the above 
requirements, and thus one can approximately present 
the tensor D in the following form: 

ro Ai(ai,a2,k) = 

Di {aiia2j + a2iaij - 2 (ai • a2) %/3 ) + 
(hkj - %/3)[i?2 (ai • aa) -t- 1?3 (ai • k)(a2 • k)] + 
£'4[(ai • k)(a2ifcj -1- a2jki) + (a2 • k)(aiifcj + ay fc^)- 

45,y (ai-k)(a2-k)/3]. (A5) 

Note that we have introduced a characteristic scale pa- 
rameter, molecular breadth ro which is used to make the 
constants -Di_4 dimensionless. 

In order to evaluate of the constants Di_4 it is conve- 
nient to convolute the tensor (|A5P with another second 
rank tensor. There are several possible second-rank ten- 
sors which can be constructed from the vectors ai , a2 
and k. We use the simple tensor aiia2j — Sij, which is 
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FIG. 7; Coefficients D1-4 for the tensor D|| as functions of 
the dipole location on the molecular long axis. The molecular 
axial ratio is set to be equal to 4. 

orthogonal to the long axes of both molecules in the case 
of ideal nematic order. It results in a convenient form of 
the convolution product: 

Z)(ai,a2,k) = ro D :(ai a2 - I) = 
+ Di [(ai • kf + (a2 • k)^] + (ai • aa)^ [Di/3 - D2/3] + 
(ai • a2)(ai • k)(a2 • k)[D2 - D3/3 + 2Di/3] + 

Z?3(ai-k)2(a2-k)2. (A6) 

Eq. (|A6|) can be expressed as a function of the molec- 
ular tilting angles 71.2 and the difference of molecular 
azimuthal angles (p — ip2 — fi'. 

£'(71,72,0) ^ do + di [P2(cos7i) + P2(cos72)] + 
d2 P2(cos7i) P2(cos72)+d3 siu^ 71 sin^72 cos2(?!)-|- 
di sin 271 sin 272 cos0 (AT) 

where the coefficients are 

do = lODi/g - 27:>2/27 + 20i:'4/27, (A8) 

di = 2L>2/9 + 47:>3/27 + 22i:»4/27, (A9) 

d2 = 2L>i/9 + 2L>2/9 + 8D3/27 + 8Di/27, (AlO) 

d3=Dj6-D2/6; (All) 

d4 = £'i/6 + i:'2/i2-i:'3/i2 + i:'4/6. (A12) 

It is possible to determine the coefficients di^^ by cal- 
culating numerically the coefficients of the spherical har- 
monic representation of D (see e.g. for details). Then 
the coefficients I?i_4 can be expressed as 

Di = -di + 5^2/4 + 5^3/3 -t- 8di/3, (A13) 

D2 = -di 5^2/4- 13^3/3 + 8^4/3, (A14) 

D3 = 9^2/4 + d3- 4di, (A15) 

1)4 = 3di/2- 3^2/4 + d3- (A16) 

In Figure[7]the coefficients -Di_4 for D|| are presented. 
They have been calculated for the molecules with ellip- 
soidal core. One can readily see that the coefficient D3 
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clearly dominates in D, and thus one may neglect other 
coefficients and write 

^TT y ~ (ai • k)(a2 • k){k,kj - %/3). (A17) 



Obviously, the approximate form (jASp changes sign un- 
der the transformations (|A4p . which means that D|| = 
— D||. As a result one can write the integrated interac- 
tion energy as 



over all orientations of the molecule '2', ?7mf(1) 
p (t/(l,2))2, which yields 



<f(1) = D, sin^ 

[2 



a X 



|Gfc(ai.k)(bi.k)-ir(ai.k)(bi.t 



, (A21) 



^^(1,2) = / C/^(l,2) d^R = 

ifj,+ -(i-)-I)^^-{fj,+ -(i-) (A18) 



and 



Substituting the dipole moments (|36p one arrives at the 
potential which depends on both long and short molec- 
ular axes. The uniaxial part which depends only on the 
long axes provides a small correction to initial uniaxial in- 
teraction potential, and we omit it here. Secondly, there 
is a part which depends on a coupling of the short axes 
of the two molecules: 



U^^{1,2) = ^i^3sin2a(ai •k)(a2 •k)x 

(bi.k)(b2-k)-i(bi.b2) 



(A19) 



Finally there exists the most important contribution 
which contains terms describing the coupling of the short 
axis of one molecule with the long axis of the other. This 
part is responsible for the induction of biaxial ordering 
by the tilt of long axes which is of primary importance 
for the description of ferroelectricity in the smectic C* 
phase. This part of the potential can be written in the 
form 

U^^(l^2) = ^^D3 sin2a(ark)(a2-k)[(bi-k)(a2-k)-H 

(ai • k)(b2 • k) - i(bi • a2) - i(ai • b2)], (A20) 

The corresponding contribution to the mean-field po- 
tential is obtained after averaging of (jA19p and (jA20p 



<f(1) 



9ro 



p Z?3 sin 2a x 



[(25 -I- l)(ai • k)(bi • k) + Gfe(2P2(cos7) -I- 1)- 

^t/(ai • k)(bi -t) - ^rsin27Cos^], (A22) 

where the order parameters Gk = ((a • k)(b • k)) and 
r = ((a • k)(b • t)) have been introduced. Apparently, 
the terms containing V (ai • k)(bi • t) and F sin 27 cos (/3 
induce the biaxial order parameter F below the tilting 
transition. 

The order parameter Gk is of minor importance, since 
it is nonzero already in the Sm^ phase. We have found 
that it is normally of the order of 0.1 and does not af- 
fect the transition significantly. The biaxial potential is 
significantly simplified if one neglects all terms contain- 
ing Gk and (a • k) (b • k) . The corresponding simplified 
mean-field potential can then be written as 



Umf = W5[cota y(a- k)(b • t)-f 

cot a F sin 27 cos ip + V2T (a • k) (b • t)] , ( A23) 



where the parameter = — ■\/2/3 pp^r^^^D^ sin^ a is 
negative and can be easily estimated for given values of 
molecular breadth rg, 2D molecular number density p, 
dipole strength /x, and dipole orientation angle a. The 
reasonable values of the constant are between 4 and 
8 (see Fig. El). 



